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The time evolution of magnons subject to a time-dependent microwave field is usually described 
within the so-called "S-theory", where kinetic equations for the distribution function are obtained 
within the time-dependent Hartree-Fock approximation. To explain the recent observation of 
"Bose-Einstein condensation of magnons" in an external microwave field [Demokritov et al, Na- 
ture 443, 430 (2006)], we extend the "S-theory" to include the Gross-Pitaevskii equation for the 
time-dependent expectation values of the magnon creation and annihilation operators. We explicitly 
solve the resulting coupled equations within a simple approximation where only a single condensed 
mode is retained. We also re-examine the usual derivation of an effective boson model from a realis- 
tic spin model for yttrium-iron garnet films and argue that in the parallel pumping geometry (where 
both the static and the time-dependent magnetic field are parallel to the macroscopic magnetization) 
the time-dependent Zeemann energy cannot give rise to magnon condensation. 
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I. INTRODUCTION 

When ordered magnets are exposed to microwave radi- 
ation of sufficiently high power, one typically observes an 
exponential growth of the population of certain groups 
of spin-wave modes during some intermediate time in- 
terval. This is an example for a general phenomenon 
which is usually referred to as parametric resonance. 
A particularly suitable system for observing parametric 
resonance are yttrium-iron garnet (YIG) crystals, be- 
cause the spin-waves in this system have a very low 
damping.^ Early microscopic theories explaining para- 
metric resonance in magnetic insulators have been de- 
veloped by SuhlP and by Schlomann and co-authorsP 
In the 1970s Zakharov, L'vov, and Starobinets^ devel- 
oped a comprehensive kinetic theory of parametric res- 
onance in magnon gases which is sometimes called "S- 
theory" . In this approach kinetic equations for the time- 
dependent distribution functions nk(t) = (a^(f)afc(t)) 
and Pk{t) = (a-k(t)ak(t)} are derived within the self- 
consistent time-dependent Hartree-Fock approximation. 
Here a& (t) and a k (t) are the annihilation and creation op- 
erators of magnons with momentum k in the Heisenberg 
picture. Subsequently the non-linear kinetic equations of 
the "S-theory" and extensions thereof have been studied 
by many authorsP^ 

Quite recently Demokritov and co-worker d 11 * 12 ! ob- 
served a new coherence effect of magnons in YIG un- 
der the influence of an external microwave field which 
they interpreted as Bose-Einstein condensation (BEC) 
of magnons at room temperature. A similar phe- 
nomenon has been observed in superfluid 3 He, where 
NMR pumping can cause the magnetization to precess 
phase-coherently.^ The emergence of this c oherent state 
can also be viewed as magnon BECPU^I Whether or 
not the experiments by Demokritov et a?! 11 * 12 ! can be 
considered to be an analogue of BEC in atomic Bose 
gases (which nowadays is routinely realized using ultra- 



cold atoms in an optical tra p) ha s been discussed con- 
troversially in the literature! 16 * 17 ! We argue below that 
the coherent state generated in these exp eriment d 1 1 * 1 2 1 
should perhaps not be called a Bose-Einstein conden- 
sate, because the condensation is not accompanied by 
spontaneous symmetry breaking in this case; instead, 
the microwave field gives rise to a term in the hamil- 
tonian which explicitly breaks the £/(l)-symmetry of the 
magnon hamiltonian. 

Unfortunately, the conventional "S-theory" is insuffi- 
cient to describe the experimental situation, because the 
coherent magnon state generated in the experiments is 
characterized by finite expectation values of the magnon 
annihilation and creation operators ak(t) and aj^i) for 
certain special values of k. In the condensed phase, the 
kinetic equations for the pair correlators nk (t) and pk {t) 
should therefore be augmented by equations of motion 
for the expectation values (ak(t)) and (aj^i)). Recall 
that in the theory of the interacting Bose gas the corre- 
sponding equation of motion for the order-parameter is 
called Gross-Pitaevskii equation}^ this equation is miss- 
ing in the conventional "S-theory" which therefore does 
not completely describe the coherent magnon state in 
the regime of strong pumping. In this work we shall out- 
line an extension of "S-theory" which includes the order 
parameter dynamics on equal footing with the kinetic 
equations for the distribution functions. Since we would 
like to clarify conceptual points rather than performing 
explicit quantitative calculations, we shall derive our ex- 
tended "S-theory" within the framework of a simple toy 
model which we motivate in the following section. 



II. TOY MODEL FOR PARAMETRIC 
RESONANCE IN YIG 

In order to understand a complex physical phe- 
nomenon, it is sometimes useful to study a simplified 
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"toy model" which still contains some essential features 
of the phenomenon of interest. For our purpose, it is suf- 
ficient to consider a single anharmonic oscillator with an 
additional time-dependent term describing the creation 
and annihilation of pairs of particles. The hamiltonian is 

H{t) = e aW ye-'^'aV + |L^ ot aa 

+ —a)a)aa. (1) 

Here a and a' are bosonic annihilation and creation op- 
erators, eo > is some energy scale, and u > is 
the interaction energy. The second and third terms 
on the right-hand side of Eq. ([I]) describe the effect 
of an external microwave-field which oscillates with fre- 
quency loq > and couples with strength 70 to the 
magnon gas. Below we shall show that this model con- 
tains the essential physics of parametric resonance and 
BEC of magnons; in particular, in the regime of strong 
pumping |7o| > |eo — wo/2| the model has a stationary 
non-equilibrium state which corresponds to the coherent 
magnon state obse rved in the experiments by Demokri- 
tov and co-workers! 11 * 12 ! 

Our toy model involves only a single boson operator 
representing the magnon at the minimum of the disper- 
sion which is expected to condense. Of course, for exper- 
imentally relevant macroscopic samples of YIG a more 
realistic model should describe infinitely many magnon 
operators labeled by crystal- momentum fc, so that the 
following bosonic "resonance hamiltonian" should give a 
better description of the experimental situation, 

H les (t) = efca^afc 

k 

k 

+ \ u ( k ^ k '^) a k+ q a k'- q a k'ak- (2) 

k,k' ,q 

If we assume that the k = boson condenses and re- 
tain only this degree of freedom on the right-hand side 
of Eq. (2), we arrive at our toy model |l]). In the the- 
ory of superfluidity a similar reduced description involv- 
ing only the order parameter is provided by the Gross- 
Pitaevskii equation.^ Of course, the minimum of the dis- 
persion in experimentally relevant samples of YIG occurs 
at certain non-zero wave- vectors ±fc*, so that it would be 
more accurate to retain the two modes a^, and a_fc >> and 
their mutual interactions in Eq. Q. Moreover, the fact 
that in the experiment ;) 11 ! 12 ! the wave- vectors of the con- 
densed magnons are different from the wave-vectors of 
the magnons which are initially generated by microwave 
pumping cannot be described within the framework of 
our toy model. Nevertheless, below we shall show that 
our simple model allows us to understand some concep- 
tual points related to the nature of the coherent state 
observed in the experiment s ! 1 1 * 1 2 1 



The bosonic resonance hamiltonian ([2]) has been the 
starting point of several theoretical investigations of 
parametric resonance in magnon gases .EHUD This model 
is believed to be a realistic model for YIG in the par- 
allel pumping geometry, where the static and the time- 
dependent components of the external magnetic fields are 
both parallel to the direction of the macroscopic magneti- 
zation. In the appendix we shall critically re-examine the 
usual derivation of Eq. ([2| from an effective spin hamil- 
tonian for YIG and show that in spin language the time- 
dependent resonance term in the second line of Eq. ^ 
involves also the combinations cos(u)ot){Sf Sf — Sf Sf] and 
$m(wot)[Sf Sf + SfSf], where S" are the components of 
the spin operators at lattice site i. Terms of this type 
cannot be related to the Zeemann energy associated with 
a time-dependent magnetic field parallel to the magne- 
tization. This is obvious for a ferromagnet with only 
exchange interactions, because in this case the magnon 
operators and a k can be identified with the Fourier 
components of the Holstein-PrimakofP^ bosons a, and 
a\ , which in turn can be related to the usual spin ladder 
operators and S~; to leading order for large spin S, 

S+ « V2S ai , Sr « V2Sa\. (3) 

Note, however, that the spin Hilbert space has only 25" + 1 
states per site, whereas the bosonic Fock space associated 
with the canonical boson operators Oj and a\ is infinite 
dimensional; the identification of magnons with canoni- 
cal bosons is therefore only approximate. For a descrip- 
tion of coherence phenomena involving large occupan- 
cies of magnon states one should therefore keep in mind 
that there is a constraint on the magnon Hilbert space. 
Assuming for simplicity that the parameter 7fc = 7 in 
Eq. ([2]) is real and independent of fc, the second term on 
the right-hand side of Eq. ^ can be written as 

k 




1 



= ^Y,{cos("ot)[S*S?-S?S?} 

i 

-sm(u;ot)[SfSf + S^Sf]}. (4) 

In spin language, the pumping term in Eq. ([2| therefore 
corresponds to a time-dependent single ion anisotropy 
whose easy axis rotates with frequency wo around the 
z-axis. Of course, the magnon operators for YIG are 
not directly related to Holstein-Primakoff bosons because 
an additional Bogoliubov transformation is necessary to 
diagonalize the quadratic part of the boson hamiltionian. 
Nevertheless, we show in the appendix that also in this 
case the pumping term in the effective boson hamiltonian 
([2]) can be related to a rotating easy axis anisotropy of 
the above type. 
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III. KINETIC EQUATIONS 

To discuss the time evolution of our toy model de- 
fined in Eq. it is convenient to remove the explicit 
time dependence from the hamiltonian H (t) by perform- 
ing a canonical transformation to the "rotating reference 
frame" , 



us explicitly give the solution of Eqs. ( 12a|12b I with ini- 
tial conditions n(0) = no and p(0) = 0. For simplicity, 
we assume in the rest of this work that 70 is real and 
positive; the case of complex 70 = |7o|e* v can be reduced 
to real 70 > by absorbing the phase factor e %Lp into a re- 
definition of the anomalous correlator, e~ %v p{t) — > p(t). 
Defining 



a = e^ ot a = U {t)alJl{t), 
a f = e-^at = U {t)(JU ] Q {t), 



(5a) 
(5b) 



where Uo(t) — e a Uoto The new operators satisfy the 
Hcisenberg equations of motion 



id t a — [a,H] , idt<v — [a\H\, 



(6) 



where the rotated hamiltonian H of our toy model does 
not depend explicitly on time, 

H = £Qa> a + —a> a> + —aa + —a> a ] aa. (7) 

Here we have introduced the shifted oscillator energy 

w 



£o = eo 



(8) 



To relate correlation functions in the original model to 
those in the rotating frame, we simply have to insert the 
appropriate phase factors. For example, in "S-theory" 
one usually considers the normal distribution function, 

n(t) = (at(t)o(t)> = (ot(t)o(t)>, (9) 

and its anomalous counter-part, 

p{t) = (a(t)a(t)) = e-^ t (a(t)a(<)> = e"^*), ( 10 ) 

where expectation values are with respect to some density 
matrix p(to) specified at time to, 



(...)= Ti-[p(t ) . 



(11) 



Throughout this work we shall mark all quantities defined 
in the rotating reference frame by a tilde. 



A. Instability of the non-interacting system 

In the non-interacting limit (u = 0) the equations of 
motion for the distribution functions n(t) and p(t) can be 
obtained trivially from the equations of motion ^ of the 
operators a(t) and c^(i) in the rotating reference frame, 

id t n[t) = 7or(*)-7oP(*). (12a) 
id t p{t) = 2e p(i)+7o[2n(t) + l]. (12b) 



a = \lt 2 o - To, 



(13) 



the solution in the weak pumping regime 70 < |eo| can 
be written as 



Rep(t) 
n + \ 
lnrp(i) 
n + \ 
n(t) + 1 



_ 1 — cos(2crf) 
-7o£o n : 



or 
sin(2ai) 



n 



-To- 



ol — cos(2at) 

l+7o j 

or 



(14a) 
(14b) 
(14c) 



In the opposite strong pumping regime 70 > |eo| the so- 
lution can be obtained by replacing a — > i/3 in the above 
expressions, where 



P = yf- 



7o 2 



Then we obtain 

Rep(t) 



n + 5 
Inrp(i) 
n Q + \ 
n(t) + \ 



= -To- 



ne 



cosh(2^) - 
-Toeo 

sinh(2/3i) 
, cosh(2/3i) 



7o" 



The behavior at the threshold value 70 = 
obtained either from Eqs. (14a 14c I for a 
Eq. ( |T6a|fl6c| ) for /3 -s- 0, 

Rep(i) 



n + 2 
Imp(i) 
n + \ 
n(t) + \ 



-2 7o ? * 2 , 



= -270*, 



n 



1 



2 7o 2 i 2 . 



(15) 

(16a) 
(16b) 

(16c) 

I ?o I can be 
■ 0, or from 

(17a) 
(17b) 
(17c) 



Physically, the exponential increase of correlations for 
To > |?o| is a consequence of the fact that in this regime 
the non-interacting part of the hamiltonian H in Eq. ^ 
is not bounded from below. This is easily seen by setting 



X + iP 



X-iP 



(18) 



These equations can be solved exactly. For |e | > | To I 
the solution is oscillatory, while in the strong pumping 
regime I70I > |?o| the solutions grow exponentially. Let 



so that 

_ 70 r~+~t , ~~1 

eoa' a + — [a 1 a 1 + aa\ 



£ ~ 70 p2 



e ^±Jl X \ (19) 
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Obviously, for 70 > |e | the non-interacting part of our 
toy model describes a harmonic oscillator with negative 
mass. The spectrum of such a quantum mechanical sys- 
tem is not bounded from below, which gives rise to the 
exponential growth of correlations discussed above. For- 
tunately, this pathology of the non-interacting limit is 
cured for any positive value of the interaction. The phys- 
ical consequences of this are most transparent if we con- 
sider the equations of motion for the expectation values 
of the creation and annihilation operators, which will be 
discussed in the following subsection. 



B. Gross-Pitaevskii equation 

The toy model hamiltonian Q in the rotating refer- 
ence frame gives rise to the following Heisenberg equation 
of motion for the annihilation operator, 



id t a = e a + 7 a^ + ua'a 2 . 



(20) 



Taking the expectation value of both sides and factorizing 
the expectation value of the interaction term as follows, 



(aV) (^}(a) 2 , 



(21) 



we obtain the Gross-Pitaevskii equation for the time- 
dependent order-parameter <p(t) = (a(t)} in the rotating 
reference frame, 



id t 4> = e o + 70^* + u\cj>\ 2 4> = 



dH c: 



(22) 



where the effective classical hamiltonian H c \ is given by 
H cX (cj>*,ct>) = e o |0| 2 + ^W 2 + cb 2 } + ^| 4 . (23) 
Writing 0= (X + iP)/y/2 we may alternatively write 



H cl (X, P) 



£ o - 70 p 2 



£0+70 ^2 



;(X 2 +P 2 



(24) 

Because the classical hamiltonian H c i(X(t), P(t)) is con- 
served along the flow defined by the Gross-Pitaevskii 
equation, the solutions of Eq. ( p2| are simply given by the 
curves of constant H c i(X(t),Pjt)) in phase space. The 
shape of H c \ and typical trajectories are shown in Fig. [l] 
Note that in the strong pumping regime 70 > |?o| the 
function H c \ (X, P) has two degenerate minima at 



X = , P = ±R = ± 



2(7o - co) 



(25) 



corresponding to stationary points (in the rotating ref- 
erence frame) of the system. Note that at these special 
points the expectation value of the annihilation operator 
is purely imaginary, 



(a) 



±- 



-P* 



7o - eo 



(26) 



(a) X/|P*| , 
1 



H ol (X,P) I u 




(b) 



H d (X,P) / u 



FIG. 1: (Color online) Graph of the classical hamiltonian 
H C \(X,P) defined in Eq. (241. The corresponding classical 
hamiltonian equations of motion are equivalent to the Gross- 
Pitaevskii equation ( 22 1 for the complex order parameter 
4>{t) = (X{t)+iP{t))JV2. The thick black lines are solutions 
of the equations of motion for different initial conditions. X 
and P are both measured in units of the momentum scale 
|P*| = \/2j7o — £q\/u. (a): eo/u = 10 and = 2; note that 
for I ?o I > 7o our classical hamiltonian H c \ (X, P) has a global 
minimum for X = P = 0. (b): eo/u = 10 and 70/w = 40; 
in the regime 70 > | e 1 our classical hamiltonian has two de- 
generate minima at (X,P) = (0,±P*), so that the graph of 
H c i(X,P) shown in (b) has some similarity to the shape of 
Napoleon's hatP 1 



The associated stationary points of the dynamical system 
( 22 1 describe a coherent magnon state where the macro- 



scopic magnetization has a rotating component perpen- 
dicular to the static magnetic field. In bosonic language, 
such a state corresponds to a coherent s tate, which is 
an eigenstate of the annihilation operator! 10 ! 22 ! Whether 
or not this state should be called a Bose-Einstein con- 
densate of magnons seems to be a semantic question. In 
our opinion this terminology is somewhat misleading, be- 
cause this coherent magnon state does not exhibit spon- 
taneous symmetry breaking which is one of the most im- 
portant properties of a Bose-Einstein condensate in in- 
teracting Bose gases. Instead, the coher ent m agnon state 
observed by Demokritov and co-worker d 11 ! 12 ! is generated 
by an external pumping field which explicitly breaks the 
[/(l)-symmetry of the magnon hamiltonian. In the static 
limit, the role of a similar symmetry breaking term on 
the Bose-Einstein condensation of magnons has recently 
been discussed by Dell'Amore, Schilling, and KramerPS 
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C. Time-dependent Hartree-Fock approximation 

Let us now take into account the leading fluctuation 
correction to the replacement (211 in the derivation of 



the Gross-Pitaevskii equation ( |22j ). To first order in u, 
fluctuations simply renormalize the bare parameters £o 



and 70 in Eq. (22 1 as follows, 



e -> e c (t) = e + 2un c (t), 
7o -> 7c(*) = 7o + wp c (t), 



(27a) 
(27b) 



where the connected correlation functions n c (i) and p c (i) 
in the rotating reference frame are defined by 



n c (t) = (Sa\t)Sa{t)), 
p c (t) = (Sa(t)5a(t)), 



(28a) 
(28b) 



with 5a(t) — a(t) — {a(t)}. Instead of the Gross-Pitaevskii 
equation ( 22 1 we now obtain for the order parameter dy- 
namics, 



id t cj> = e c (t)<f> + 7 c (t)^* + u\cj>\ 2 (j). 



(29) 



Note that this generalized Gross-Pitaevskii equation de- 
pends on the connected correlation functions n c (t) and 
p c {t), which we calculate in self-consistent Hartree-Fock 
approximation. The resulting equations of motion can be 
obtained from the corresponding non-interacting kinetic 
equations ( 12a|12b| by substituting 



eo->e(t) = i Q + 2u[n c {t) + \<P(t)\ 2 } 
lo^ l(t) = jo + u[p c (t) + cj) 2 (t)}. 



(30a) 
(30b) 



The kinetic equations for the connected distribution func- 
tions are therefore 



id t n c (t) = -y(t)p*(t) - 
idtpc(t) = 2e(t)p c (t) 



l*(t)Pc(t), 

- 7 (t)[2n c (t) 



!]■ 



(31a) 
(31b) 



For cp — these equations reduce to the kinetic equations 
obtained within "S-theory" ™ The numerical solution of 
Eqs. (|9j [31a[ |31bJ) for n c (0) = n , p c (0) = 0, and in- 
finitesimal lm<p(t) > is shown in Fig. [2] Obviously, 
for sufficiently strong pumping an infinitesimal initial 
value of 4>(0) builds up to a finite oscillation. Moreover, 
the connected correlation functions n c (t) and p c {t) re- 
main always bounded, in contrast to the exponentially 
growing correlations in the non-interacting limit given in 
Eqs. (16a 16c). Note also that the time evolution of the 



connected correlation functions appears to be rather ir- 
regular as soon as the order-parameter has built up to 
a finite value. In the conventional "S-theory" the quan- 
tities n c and p c are periodic (Fig. [2j;), while including 
the order parameter dynamics disturbes this strict peri- 
odicity (Fig. [2]d) . This feature is still missing within the 
usual "S-theory" in the strong pumping regime. 
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FIG. 2: (Color online) Numerical solution of the coupled ki- 
netic equations (29 31a 31b I with initial conditions n c (0) = 
no — 1, p c (0) — 0, and 0(0) = 0.05i. The character- 
istic energy scales a and /3 are defined in Eqs. (13 15 1. 

(a) : io/u = 500 and 70 /it = 200. Recall that in the absence 
of interactions there is no instability as long as |?o| > 70. 

(b) : eo/u = 500 and 70/w = 5000. In this regime there would 
be an instability in the non-interacting limit, but in the in- 
teracting system all correlations remain finite, (c): Same pa- 
rameters as (b) but without finite expectation values, like in 
the conventional "S-theory". 



IV. SUMMARY AND CONCLUSIONS 

Let us briefly summarize the two main results of this 
work: 

First of all, we have shown that a complete theoreti- 
cal description of the coherent magnon state emerging in 
YIG for sufficiently strong microwave pumping requires 
an extension of the usual "S-theory" which includes the 
Gross-Pitaevskii type of equation for the expectation val- 
ues of the magnon operators. Within a simple toy model 
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consisting only of a single magnon mode we have shown 
how to construct such an extension. The explicit solution 
of the resulting kinetic equations shows that the order 
parameter dynamics strongly influences the distribution 
functions. 

Our second main result is the observation that in 
spin-language the usual bosonic resonance hamiltonian 
([2]) corresponds to a time-dependent rotating easy axis 
anisotropy whose axis is perpendicular to the direction of 
the external field. If this anisotropy is sufficiently strong, 
it gives rise to a forced oscillation of the macroscopic 
magnetization around the direction of the static exter- 
nal field. Although this phenomenon can be described in 
terms of a coherent magnon state, it should not be called 
a Bose-Einstein condensate, because the emergence of 
this state is not associated with any kind of spontaneous 
symmetry breaking. 

In future work, we shall further extend our approach 
in two directions: on the one hand, a realistic model for 
YIG involves a quasi-continuum of magnon modes, which 
can condense at finite wave- vectors ±fe*. For a more 
realistic quantitative description of the experiments, we 
should therefore generalize our extended "S-theory" to 
include all magnon modes relevant to the experiments 
on YIG. This would also allow us to distinguish between 
the "primary magnons" created by the external pump- 
ing, and the "condensing magnons" with wave-vectors at 
the minima of the dispersion. The second direction for 
improving our approach is to include correlation effects 
beyond the self-consistent Hartree-Fock approximation 
into the kinetic equations. For example, to second order 
in u the kinetic equations will contain relaxation terms 
which will damp the oscillatory time dependence found 
at the Hartree-Fock level. Work in both directions is in 
progress. 



model™ 

H Y i G {t) = 



ij ap 



-[h a + hiCos,{bjQt)]^Sf, 



(Al) 



where a, j3 — x,y, z label the three spin components, and 
the exchange couplings JJy = J(r^ — rj) are only finite if 
the lattice sites rj and tv, are nearest neighbors on a cubic 
lattice with lattice spacing a ~ 12.376A. The value of 
the nearest neighbor exchange is J « 1.29K. The dipolar 
tensor D"? — D a P (r.; — rj ) is explicitly 



D 



a0 



(!-<%)] 



Of." f/ 3 



gap 



(A2) 



where r,j = r ; — rj and fij = / \ | . If we arbitrar- 
ily set the magnetic moment /1 = 2 /is — eh/(mc), then 
we should work with an effective spin S « 14.2, as dis- 
cussed in Ref. [21]. Here ho and hi are the amplitudes 
of the static and oscillating magnetic field (multiplied by 
/i). We assume that ho > \hi\ and that both the static 
and the oscillating magnetic field point into the direction 
of the macroscopic magnetization which we call the z- 
axis. At this point one might already wonder how in this 
parallel pumping geometry one can possibly arrive at a 
bosonic resonance hamiltonian of the form ^ , which ac- 
cording to Eq. Q can be related to some rotating easy 
axis anisotropy. In fact, we shall show shortly that the 



spin hamiltonian (Al) with parallel pumping cannot be 
reduced to the bosonic resonance hamiltonian ([2]). 



To bosonize the hamiltonian ( Al ) we express the spin 



operators in terms of boson operators bi and b\ by means 
Holstein-Primakoff transformation's 
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APPENDIX: PARALLEL PUMPING OF 
MAGNONS IN YIG 

It is generally accepted that the magnetic proper- 
ties of YIG in the parallel pumping geometry can be 
modelled by the following time-dependent quantum spin 



= 1 - g&i = ( A3a ) 



S-b\h. 



(A3b) 



As usual, the square roots are then expanded in powers 
of 1/5, resulting in a hamiltonian of the form 



H YlG (t) = H (t) + H 2 (t) + H h 



(A4) 



where H (t) is a time-dependent constant, H 2 (t) is 
quadratic in the boson operators, and the time- 
independent interaction Hi nt involves three and more bo- 
son operators. After Fourier transformation to momen- 
tum space the quadratic part of the hamiltonian can be 
written as 



A k b{b k + ^fblbl k 



B. 



Ht Lt 



cos(w £) b\b k 



B 



rb- k b k 



(A5) 
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where 



with 



i A 



Bk — B-f, — y ^2 l e lhr,J Bij, 



Aij — SijflQ + S(Sij S Ji n — Jjj) 



5 



rjza; _i_ T)VV 



Si 
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(A6a) 
(A6b) 



(A7a) 
(A7b) 



Finally, we use a Bogoliubov transformation to diagonal- 
ize the time-independent part of H2 (t) , 



Uk ~Vk 

-v% u k 



(A8) 



showed that the neglected terms can qualitatively change 
the results obtained in resonance approximation. Here 
we would like to point out that the approximations lead- 
ing to Eq. (A13l amount to an essential modification of 
the original spin hamiltonian. To see this, let us for the 
moment accept the validity of these approximations, thus 
replacing Eq. ( All[ ) by the non-interacting part of the 
resonant hamiltonian ([2]), 



^ £ka k a k 

k 



in 



l k e - iuot a{al k + -fie**** a_ k a k 



, (A14) 



where we have dropped the constant terms. Using now 



the inverse of the Bogoliubov transformation ( A8 ) to re- 



express the magnon operators in Eq. (A14| in terms of 



Holstein-Primakoff bosons and assuming for simplicity 
that 7fc is real, the second term in Eq. (A14| can be 
written as 



where 



A k + £fc 



2e fe 



Vk 



Bk 
\B k \ 



Ak — £fc 



2e fc 



and 



(A9) 



(A10) 



After this transformation the hamiltonian reads^ 



tka\ak 



£fc — Ak 



+hi cos(u> t) £ 



.4 



(k 
,t„t 



a h a,k + 



A k 



2e k 



+ £ [ 7fc cos ( w oi)aLaI fc + Ik cos(o; i)a_ fc a fe 
k 



where 



Ik 



h\Bk 
' 2e fe ' 



(All) 



(A12) 



To obtain the quadratic part of the resonance hamilto- 
nian (J2J) from Eq. (All) two additional approximations 
are necessary: the second line on Eq. (All I involving the 



combination cos(uiot)Aka k ak has to be dropped, while in 
the last line one should substitute 

7fc cos( Wo i) -> ^e-^ ot , 7 feCOsM) -> ^e^*.) 

(A13) 

Apparently this approximation has been accepted for 
many decades in the literaturel^ES] However, a thorough 
study of the non-resonant terms neglected in this ap- 
proximation has been performed by Zvyagin el al. pi who 



1 ^(jkAk cog ^ o ^ 
k ^ £k 



b l b -k + b -kbk 



+ijk sin(w t) 

ElkBk 1 ,x 
costal) 
e k 



b k b -k - b -k b k 



blbk 



(A15) 



Only the last term on the right-hand side has the form of 
the boson representation of the Zeemann term associated 
with an external pumping field parallel to the magneti- 
zation, while the first two terms can be identified with 
the boson representation of spin anisotropics associated 
with a rotating easy axis perpendicular to the z-axis, 
see Eq. Q. We thus conclude that the time-dependent 
part of the resonant hamiltonian ^ does not represent 
the time-dependent Zeemann energy associated with a 
harmonically oscillating magnetic field in the direction 
of the magnetization. Instead, the time-dependent off- 
diagonal pumping terms arise from a rotating easy axis 
anisotropy perpendicular to the magnetization. The mi- 
croscopic origin of such a term is not clear to us; pos- 
sibly the time-dependent electric field associated with 
the harmonically varying magnetic field parallel to the 
magnetization can indirectly induce such a term in the 
spin hamiltonian, similar to the second order interaction 
hamiltonian in the theo ry of t wo-magnon Raman scatter- 
ing in antiferromagnets! 25 * 26 ! Moreover, in real materials 
crystallographic or shape anisotropies can give rise to fur- 
ther contributions to the effective spin hamiltonian which 
after Holstein-Primakoff transformation might have the 
same form as the terms in Eq. (A15). 
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